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we study the renormalization group flow of the gauge coupling constant. Our calculation gives 
consistent result with the perturbative prediction of the running coupling in the weak coupling 
region. We investigate a possible signal for conformal behavior in the strong coupling region. 
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1. Introduction 

There has been a lot of interest in the non-supersymmetric gauge theories with non-trivial 
infrared fixed point (IRFP) [|l|, Such a theory could play an important role in the phenomenology 
of the physics beyond the standard model. Since the existence of the IRFP gives slow running of 
the gauge coupling, it may be a possible candidate of walking technicolor model ^, ||, |6|, 
Within the framework of the walking technicolor scenario, the large anomalous dimensions for 
the bilinear operators could give a possible way to resolve the tension between the realistic quark 
masses and the suppressed flavor changing neutral current. 

SU(N) gauge theories with many flavors have been considered as traditional candidates for 
the walking technicolor scenario. In fact, the two-loop perturbative j8 function predicts the exis- 
tence of the IRFP within a certain region of the number of fermions which is so called conformal 
window. Near the conformal window, the anomalous dimension is expected to be enhanced from 
the Schwinger-Dyson analysis. In addition, in such a theory the contribution to the electroweak S 
parameter might be suppressed in contrast to QCD like theories. These features are favorable from 
the view point of model construction of the realistic dynamical electroweak breaking. Furthermore, 
once a theory has an IRFP, it does not break the chiral symmetry and should be inside the conformal 
window. It is helpful to understand the phase structure of the gauge theory, i.e. whether a model is 
in the chiral symmetry broken phase or not. The identification of the conformal gauge theory itself 
is interesting beyond phenomenological issues. 

In the lattice calculation, whether a certain model has an IRFP or not can be answered from 
first principle. A simple way to identify the conformality is an analysis of the renormalization 
group flow of the gauge coupling constant. Recently non-perturbative lattice calculations of the 
running coupling constant are widely investigated. There are several methods in analysis, e.g. the 
Schroedinger functional (SF) running coupling scheme [||, ^] and twisted Polyakov loop method. 
There has been lots of studies of the running coupling constant e.g. rif = 8, 12 and 16 flavor models 



in SU(3) fllOl hin, two flavors of adjoint fermions in SU(2) gauge group [12], and for a recent study 



of the SU(2) with fundamental six-flavor model, see Ref. [|13|]. 

In this report, we study the infrared behavior in Uf = S flavor SU(2) model. We employ twisted 
boundary conditions in which the methods for calculation of the renormalized coupling can be de- 
fined. This twisted boundary condition does not need the additional boundary terms in contrast to 
SF boundary so the measured values have no 0{a) contribution. Therefore we could obtain sys- 
tematically reliable results. Using the finite volume scaling method, we obtain the renormalization 
group flow of the gauge coupling constant and study its low-energy behavior. In order to see the 
non-perturbative features of the theory, we investigate the possible signals for the IRFP by using 
the growth rate of the running coupling constant. This quantity is easy to compare with the one 
obtained from the perturbative calculations. 

This report is organized as follows. In section 2, we introduce two methods for calculation of 
the running coupling constant in twisted boundary conditions; one is obtained from the Polyakov 
loop correlation functions (TPL scheme) and another is from the Wilson loop. In section 3, we 
show some numerical results on the running coupling constant in 8-flavor SU(2) gauge model. In 
section 4, we show our preliminary study of the growth rate of the running coupling constant for 
each of the renormalized coupling. Finally in section 5, we summarize our results. 
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2. The running coupling in twisted boundary condition 

In this section, we present the definition of the twisted boundary conditions for both of link 



and fermion variables in SU(2) gauge theory. Here we consider twisted boundary conditions [ |14| ] 
for the link variables in x and y directions on the lattice; 

U^ix + vL/a)=ayU^Q.l, {v=x,y) (2.1) 

where the twist matrices Q.^ are introduced as follows, 

a^^y = -Q.yQ.x, O.^^l = 1, = -1, TrH^ = 0. (2.2) 

The consistency of the twisted boundary condition for fermions should accompany with the 
additional 'smell' flavor degrees of freedom. We take the following boundary conditions 

rai^ + vL/a) = e'^/'afwp i^Dpa ' (2-3) 

for V =x,y directions. Here the indices of a, j8 and a,b represent the smell degree A'^^ and color Nc- 
The number of smell should be multiples of the number of color. We take Ns = Nc = 2. 



Then the TPL scheme coupling [ ]15[ ] is defined as 



2 _ 1 {LyMy,z,L/2a)PM0,0V) 
^^'"^ k{^,,.P,ix,y,L/2a)PM0,0y)' ' 

where P^,? are the Polyakov loop in the twisted and untwisted directions defined in a gauge invariant 
way. For instance is given by 

P,{y,z,t) = Tr |^nUx(x = j,y,z,t)a,e''^y/L^ . (2.5) 

The coefficient k is obtained by analytic calculation which gives k = 0.0172 • ■ • ^ 

2.1 Method using Wilson loop 

We also introduce the method for the coupling constant by using Wilson loop [16]. We first 
use the Creutz ratio in a finite volume iLjdf' as 

where W{R/a,T /a) means the rectangle Wilson loop. Each ofR/a and T /a means the length of 
the Wilson loop for each different direction. To fix the scheme we consider a measurement of the 
Wilson loop for untwisted-untwisted directions. Taking different directions such as twisted-twisted 
direction corresponds to the different scheme of defining of the renormalized couplings. Then this 
quantity is proportional to the gauge coupling at leading order. By multiplying the tree level 
matching factor \/k, renormalized gauge coupling constant can be defined by 

gl{L/a,R/L+\/2) = x{R/a + l/2,L/a)/k. (2.7) 



' This value differs from its original value in Ref . | p^ ] because of inclusion of the trace factor of the denominator. 
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Here is a function of the ratio {R/L+\/2) = r, L/a and physical volume L(= where ji 
is the energy scale. In order to obtain the renormalized coupling g|, in a finite volume, we fix the 
value of r to an appropriate value. In our case, we take r = 0.25. 

Since xir^Lja) itself has a systematic artifact from discretization [16, 11], it is better to im- 
prove g^ij) by using discretized value of the leading factor k{r) instead of its continuum value 
which is calculated in a lattice perturbation [|n|]. While both the kir) and xir) contain large lat- 
tice artifacts, some of them may be canceled each other and the ratio x{r)lk{r) could behave as a 
smooth continuous function of r. 



3. Simulation detail 



We carry out the measurements of the coupling from TPL and Wilson loop methods. Our sim- 
ulations are performed on L = 6, 8, 10, 12, 14, 16, 18 lattices with plaquette gauge action with one 
massless staggered fermion with twisted boundary condition which corresponds 8-flavor fermions. 
The gauge configurations are generated by the Hybrid Monte Carlo algorithm. We take several 
values of j8 = 1.3 ~ 15 for each of lattices. We accumulate about 2 ~ 5 x 10"* trajectories for 
L = 6 ~ 14 and about 5,000 ~ 20,000 trajectories for L = 16 ~ 18. Then we take every trajectory 
for the measurement of the Wilson loop and Polyakov loop. The error estimation is made by a 
standard Jackknife analysis with O(IO^) bin size. 

As defined in the previous section, the TPL gauge coupling is obtained directly from the con- 
figurations as long as Lja is even. While in the case of L = 6, 10, 14, 18, we can obtain directly, 
in the case of L = 8, 12, 16, we need to interpolate the Creutz ratio as a function of r. Since k(r) 
becomes zero at r = 1/2, g^(r) blow up at r = 1/2 and we find that the function ^{f) is reason- 
ably fitted by using the pole type functions as ^(f) ~ 7^775 where f(r) is a polynomial function 
with f(r = 1/2) 7^ 0. Hence we obtain the renormahzed gauge couplings from Wilson loop method 
even in small volumes ^. 

Thus we obtain both the results of the renormalized couplings defined by two methods. The 
results are shown in Fig.[T[ As shown here, both show similar j3 dependence and they are close to 
each other in the weak coupling region. We find the relative error of the Wilson loop coupling is 
about 1/5 ~ 1/10 of that of TPL coupling even in small jS and large L/a regions. As it is important 
to know the infrared region of this model, now we concentrate on the analysis of the coupling from 
Wilson loop. 



4. Result 



In order to study the renormalization flow of the running couplings, we use the finite scaling 
method. First we define the step scaling function for the gauge coupHng as 

l.{u,s,a/L) = g^{l5,sL/a)\^2(^pLf^)^^, (4-1) 

where parameter s is the scaling parameter. In our analysis, we set s = 3/2. We obtain four data 
of g^{l5,L/a) at a fixed u for coarse lattices such that L/a = 6, 8, 10, 12 by tuning j8. In the case of 

^For detailed analysis, see 
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Figure 1: The renormalized gauge couplings for each j3 and L. In the left panel, the data show the TPL 
gauge couplings, in the right panel, the data show the Wilson loop couplings. 



Wilson loop, it is easy to interpolate to the arbitrary value of jS by using the following fit function 
as 

Our fitting procedure is as follows. First we carry out the fit of j8 dependence of gauge cou- 
plings at each of fixed lattices with the free parameters jSq and c;. We use typically 6 ~ 8 free pa- 
rameters to fit the results. Fit range of jS is taken as j3 = 1.5 15. Thus we interpolate g^{p,L/a) 
at arbitrary value of j8 for each of lattices. In order to obtain the g^{^,L/a) for unmeasured lattices 
of L/a = 9 and 15, we also fit the volume dependence of g^{^,a/L) at fixed j8 by the following 
polynomial functions 

gl.{p,L/a) =do + diL/a + d2{L/af, (4.3) 

where the coefficients dj are free parameters. In these fits we use four data points of L/a = 
6,8, 10, 12 for obtaining the value at L/a = 9 and L/a = 12, 14, 16, 18 for L/a = 15, respectively. 
Finally using four data points of Jl{u,3/2,a/L) we take continuum extrapolation and obtain con- 
tinuum step scahng function a{u) defined by a{u) = lima/L^o'^{u,s,a/L). We try to carry out the 
continuum extrapolations by two ways; one is a linear extrapolation in terms of (a/L)^ by using 
three data points ofa/L closing to the continuum limit, i.e. a/L=l/12,l/15 and 1/18 and another 
is a quadratic extrapolation in terms of {a/LY by using four data points. 

Here we show some results of the continuum extrapolation of £(m). Figure^ shows the contin- 
uum extrapolations of the 'L{u,s,a/L) with u = 0.8 (left panel) and u = 4.5 (right panel). Figure || 
shows the result of the o{u)/u. In the perturbative region the results at the continuum limit of linear 
and quadratic extrapolations o{u = 0.8) gives consistent result with perturbative 2-loop continuum 
prediction. On the other hand, in the strong coupling region. We find the result of o{u)/u are larger 
than that of two-loop continuum perturbation. Although the two results are consistent with each 
other within 1 a level and also consistent with the IRFP above g|, ~ 6, because the growth rate is 
consistent with dashed line in Fig. ^, we can not yet conclude clearly that there is an IRFP due to 
its large statistical uncertainty. In order to determine the renormalization flow of o{u), the more 
statistics are needed especially in large lattices such as a/L = 1/16 and 1/18. 
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Figure 2: The data shows the resuhs of a{u) for each of lattices L/a ~ 9, 12, 15 and 18. The continuum 
extrapolation results for linear and quadratic fits are shown by blue and red curves, respectively. The solid 
line represent the initial values of m (m = 0.8 in left panel, u = 4.5 in right panel). The dotted line represent 
the two-loop results of the a{u). 
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Figure 3: a{u)/u obtained by linear (blue) and quadratic (red) continuum extrapolation. The horizontal 
axis shows the initial value of m = g^. The vertical axis show the result of a{u)/u. For the reference we 
also show the two-loop perturbative results of the a{u)/u by solid curve. Dashed line represents a{u) = u. 



5. Summary 

In this proceedings we reported on the infrared behavior of the SU(2) gauge theory with 
fermions of 8-flavor in the fundamental representation. We measured the renormalized gauge cou- 
plings defined in twisted boundary conditions. This twisted boundary condition does not need the 
additional boundary terms so the measured values have no 0{a) contribution. Using the finite vol- 
ume scaling method, we obtained the renormalization group flow of the gauge coupling constant. 
We found that in the weak coupling region the running of gauge couplings gives consistent result 
with perturbative two-loop prediction. However, in the strong coupling we can not have conclusive 
results of the existence of the IRFP due to the large uncertainty. 

The accuracy of our results could be improved by accumulation of more statistics. This may 
help to stabilize the continuum extrapolations of the running coupling and enable us to study the 
consistency of the two results obtained by two different methods which gives non-trivial check 
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of the existence of the IRFP. We are also studying other physical quantities such as Wilson loop 
itself and fermion bilinear operators. They could also reveal some properties of the strong coupling 



gauge theories. These works are still ongoing []17[]. 
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